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We study theoretically charge transport properties of normal metal (N) / ferromagnet insulator 
(FI) / superconductor (S) junction and S/FI/S junction formed on the surface of three-dimensional 
topological insulator (TI), where chiral Majorana mode (CMM) exists at FI/S interface. We find 
that CMM generated in N/FI/S and S/FI/S junctions are very sensitively controlled by the direction 
of the magnetization m in FI region. Especially, the current-phase relation of Josephson current 
in S/FI/S junctions has a phase shift neither nor tt, which can be tuned continuously by the 
component of m perpendicular to the interface. 
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The class of time-reversal (TR) symmetric insulators 
with the nontrivial topological properties has been pro- 
posed theoretically and discovered experimentally [H, 
[3|. The hallmark of this insulator, i.e., topological insu- 
lator (TI), is the edge or surface channels, while the bulk 
states are gapped and inert for the low energy phenom- 
ena. In the two-dimensional (2D) case, there appears 
the helical modes at the edge of the sample, i.e., the pair 
of one-dimensional modes connected by the time-reversal 
symmetry and propagating in the opposite directions for 
opposite pseudo-spins (2[ . This is analogous to the chiral 
edge modes in the quantum Hall systems [4]. Therefore, 
the TI in 2D case can be regarded as the two copies of the 
quantum Hall systems for up and down pseudo-spins [5|] , 
and is often called the quantum spin Hall system. In 
three-dimensional (3D) case, there are two class of TI, 
i.e., weak TI (WTI) and strong TI (STI), corresponding 
to the even (WTI) and odd (STI) number of the chiral 
Dirac fermions on the surface [6] . Since the even number 
of Dirac fermions can be paired to open the gap, those 
in WTI is fragile against the disorder and/or the inter- 
actions while they are robust in STI. The 2D quantum 
spin Hall system is adiabatically continued to the WTI 
when the weak interlayer coupling is tuned, while the 
STI has no analogue to the quantum Hall system, and is 
a genuine new state of matter. 

The 2D chiral Dirac fermion on the surface of STI is 
protected by the bulk gap and its topological property. 
Therefore, it offers an interesting system to look for the 
2D superconductivity with the Cooper pairs mediated 
by the bosonic excitations, e.g., phonons and excitons, 
in the STI. Fu and Kane studied the superconductivity 
induced by the proximity effect to the surface of STI Q . 
Considering the interface between the ferromagnetic in- 
sulator (FI) and conventional superconductor (S), they 
predicted the appearance of the chiral Majorana state as 
an Andreev bound state [7]. We call this chiral Majorana 
mode (CMM), which has a dispersion along the interface 
while it is confined along the direction perpendicular to 
the interface. Detecting the Majorana fermions in terms 



of the interferometory has been proposed also [8|, |9| . 

The presence of CMM is predicted in the Px-\-Wy chiral 
superconductors [13], e.g., Sr2Ru04. However the con- 
trol of the chiral domains and manipulation of the edge 
channels are experimentally difficult [ll| . In addition, the 
Px + Wy superconductivity is very fragile against the dis- 
order. Also, the CMM in the ^He and cold atoms have 
been studied theoretically 0, but they are neutral sys- 
tems and charge transport is missing there. Therefore, 
the present system offers unique opportunity to study the 
quantum charge transport specific to CMM and its con- 
trol, which is more promising to be realized experimen- 
tally. However, its theoretical studies have been limited 
focusing on the detection of the Majorana fermion itself 

mm- 

In this paper, we study the manipulation of the 
quantum transport properties associated with Majorana 
fermions at the interface of the superconductor(S) and 
ferromagnetic insulator (FI) generated on the surface on 
TI. Hereafter, since we concentrate on the STI, we denote 
TI instead of STI for simplicity. We show that the direc- 
tion of the magnetization m can be used to control the 
Andreev reflection and Josephson current via the CMM 
generated in N/FI/S and S/FI/S junctions, offering a 
unique method for superconducting spintronics. 




FIG. 1: Schematic illustration of the junction, (a) Normal 
metal (N) / Ferromagnetic insulator (FI) / Superconductor 
(S) junction and (b) S/FI/S junction formed on the surface 
of 3D topological insulator (TI). The current is flowing on the 
surface of TI. 
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We consider FI/S and S/FI/S structures formed on the 
surface of 3D topological insulators as shown in Fig.l. 
We concentrate on the situation that TI below FI be- 
comes ferromagnetic insulator due to the exchange cou- 
pling. Since surface state of TI is metallic, we can regard 
the configuration shown in Fig. 1(a) as normal metal 
(N)/ FI /S junction with N{x < 0), FI(0 < x < d), smd 
S{x > d). We also consider S/FI/S junction as shown in 
Fig. 1(b) with S(x < 0), FI(0 <x <d), and S(x > d). 
The interfaces N(S)/FI and FI/S locate at x = and 
X = d^ respectively. We assume that the chemical poten- 
tial in N and S are equal to each other. 

The Hamiltonian of the surface state on TI is given by 



( H(k) + 
V -A* 



M 



A 



■M* 



(1) 



with H{k) = vpio'xkx + ^yky) — fi[Q{—x) + Q{x — d)], 
and M = m ■ &Q{d — x)Q{x) with m ■ & = rrixCFx + 
^y^y + ^z^z- Here, /i, (T, vf^ tu denote chemical 
potential, Pauli matrices, velocity, and magnetization 
(times the exchange coupling constant which we assume 
to be 1), respectively [7|. Note that m enters as an ef- 
fective vector potential A of the electromagnetic field. 
The pairing symmetry of superconductor is assumed to 
be s-wave and A is given as A = i&y AQ{x — d) and 
A = iGy[/\Q{x - d) + Ae(-x)exp(i(/p)] for N/TI/S 
and S/TI/S junctions, respectively, where denotes the 
macroscopic phase difference between left and right su- 
perconductor. In general, the magnitude of A is smaller 
than that of the bulk energy gap of superconductor de- 
posited on TI due to the nonideal S/TI interface 

First, let us consider N/FI/S junction (Fig. 
1(a)). A wave function of an electron injected 
from N with an injection angle 6 is given as 
= e-x.-p{ikyy)[^ nIx)Q{-x) + ^ Fi{x)Q(x)Q{d - 
x) + ^^{x)Q{x — d)]^ where ky = kp sinO is the momen- 
tum parallel to the interface with vrkp = /i. ^Ar(x), 
^F/(^) and (x) are given by 

^Ar(^) = [(^m + ^i^/ir) exp(z/c^x)+6^er ex.p{-ikxx)]{2) 

^F/(^)= [^el exp( — ^x) + ^e2 exp(/^*x) 

exp(^x) + '^h2 exp(-/^*x)] (3) 
(x) = [^et expiikxx) + "^ht expi-ikxx)] (4) 



with kx 



kyj K,X 



Kx + irrix/vF^ and 



Kx = Y^m^ + {vrky + rriyY /vp- The four component 
wave functions ipin^ iphr and iper are given by ^tpin = 
(l,exp(z^),0,0), ^^hr = (0,0,l,-exp(z^)), ^^er = 
(1, — exp(— i^), 0, 0), with exp(z6>) = {kx^iky) /kp- Other 
wave functions V^ei, V^e2, V^^i, V^/i2, V^et, and ^jjht are ob- 
tained by solving eq. ([T]) in a similar way. 

The coefficients of Andreev reffection a and normal 
reffection h are obtained by imposing the boundary con- 
dition ^Ar(x = 0) = ^F/(^ 0), and ^fi{x = d) = 
'^g{x = d). Then, the angle resolved tunneling conduc- 
tance for injection angle is obtained by standard way 



as cr5'(^) = 1+ I a p — I 6 p. The normalized angle aver- 
aged tunneling conductance a by its value in the normal 
state is given by [l5[ 



(J = 



r_i%<js{9) cos ede 



r'r/2 

J_^/2 COS 
CT^[1+CT^ Irp-(l-CT^) in^] 

I i + (i-(Tiv)exp(n)r2 12 



(5) 
(6) 



with exp(z7) = [rrizQ^o^O + z(/isin^ + my)]/[mzQ^o^O — 
i{lisme + my)] andF = A/{E^VE'^ - A^). aN denotes 
the transparency of the junction in the normal state given 
by (Tn = I / {cosh^ {tZxd) + tan^ 6> sinh^(/^a^(i)/c^//^^). For 
(Tat ^ 0, the denominator of as{0) becomes zero at £^ = 



A(/i sin + my)sgn{mz) 
(/i sin + rriyy + cos^ 0ml 



(7) 



This condition coincides with the formation of CMM at 
the FI/S interface with semi-infinite. CMM and An- 
dreev refiection are strongly related to each other and 
I a 1=1 is satisfied at £^ = £^5 independent of ctn. It 
is noted that the sign of £^5 is changed by reversing 
the direction of magnetization rriz. £^5, crs{^) and a 
are independent of rrix. As a special limit of rriz = /J^ 
and rriy = 0, this formula includes the case of chiral 
Px + isgn{—mz)py-wdive superconductor where £^5 is re- 
duced to be El) = —AsmOsgn{mz) [l^ although the 
given pair potential is full gap 5- wave. It is remarkable 
that the sign of corresponds to the chirality of the 
CMM, which can be controlled by the direction of mag- 
netization of FI. Here, we focus on how the above CMM 
infiuences the bias voltage V dependence of conductance 
in N/FI/S junctions with E = eV. As shown in Fig. 2, 
the resulting a has a zero bias conductance peak (ZBCP) 
originating from the peak of as{0) dit E = E^. As seen 
from the left panel of Fig. 2, the slope of the curve of 
Eb around E^ = {0 = 0) becomes gradual with the 
decrease of /i/m^. Then, the contribution around ^ = 
becomes significant in the integral of numerator in eq. 
([5]) and the resulting height of the ZBCP is enhanced 
with the decrease of the magnitude of /i/m^ as shown in 
the right panel. 

The presence of niy changes £^5 drastically since 
Eb{0) = —Eb{—0) does not hold any more. As shown 
in the left panel of Fig. 3, in the presence of m^y, E^ does 
not become zero any more for 6> = 0. For my/rriz = 0.4, 
the domain of with negative E^ is larger than that of 
positive one. The resulting a has a peak at negative bias 
voltage as seen from curve h in the right panel of Fig. 3. 
On the other hand, for rriy/mz = —0.4, the domain of 
with positive E^ is larger than that of negative one as 
seen from curve c in the left panel. Then the resulting 
(J has a peak at positive bias voltage as seen from curve 
c in the right panel of Fig. 3. As shown above, CMM, 
Andreev refiection and resulting crsiO) are critically con- 
trolled by m. The present CMM is significantly different 
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CMM Normalized Conductance 




FIG. 2: (Color online) Left Panel: Chiral Majorana mode 
(CMM) energy level Eb as a function of the incident angle 6. 
Right Panel: Normalized tunneling conductance a in N/FI/S 
junctions. rUzd/vF — 1 and rriy/mz — 0. a: fi/mz — 1, b: 
Ii/ttIz — 2 and c: ii/ruz — 0.5. 

from that in the non-centrosymmetric superconductor, 
where CMM appears as helical edge modes [17]. In the 
present case, one of the spin-split bands is missing com- 
pared with the non-centrosymmetric superconductors. 



CMM Normalized Conductance 




FIG. 3: (Color onUne) Left Panel: CMM energy level Eh. 
Right Panel: a in N/FI/S junctions. rrizd/vF — 1 and 
Ii/ttIz — 1. a: my/mz — 0, b: TJiy/mz — 0.4 and c: 
niy/mz = —0.4. 



of El) [eq. (|9])] is very general including several preexisting 
cases. As a special limit with £^5 = A, and rrix = 0, E^, 

is reduced to Eb = ±^cos2((/p/2) + (1 - (Jn) sin^((/p/2)A, 
which corresponds to the Andreev bound state in conven- 
tional S/non- magnetic insulator (NI)/S junction [19]. If 
we choose £^5 = 0, and rrix = 0, the bound state formula 
in dxy{px)-^^'y^ /NI/(ia;2yfe)~wave junction is reproduced 
p^oj . To understand the and (p dependence of ^5 {Ej) 
more intuitively, we plot (absolute value of ^5) in Fig. 
4. Ej is an even function of for any mx- For rrix = 0, 
Ej is symmetric with respect to = (Fig. 4(a)). How- 
ever, for rrix 7^ 0, the resulting Ej is no more symmetric 
around (f = as shown in Figs. 4(b) and 4(c). This fact 
seriously influences the Josephson current / given by 

Sin(^ - 2,5) de ^AHanHEjmBT)a^cos0 

bIRn = 5-7^ 

J_l/2dO<JN COS 6* 

(10) 

with resistance in the normal state Rn- 



(a) CMM (b)CMM 




Next, we focus on the Josephson current in S/FI/S 
junction (Fig. 1(b)). Since the magnitude of the pair 
potential of the left and right superconductors are equal 
to each other, it is possible to evaluate the Josephson 
current by using the CMMs formed in S/FI/S junctions 
[18^. The CMMs can be obtained as in the case of N/FI/S 
junction with wave functions of Eq. ([T]) under a proper 
boundary condition. The resulting Josephson current / 
can be obtained as 

9p E E 
I=j^, F=-..TElog[2cosh(^)],(8) 

k 

Ej = ^(JN cos2((^/2 - 5) + (1 - a^)(^5/A)2A (9) 

with 8 = rrixd/vF and temperature T. The appearance 
of S in the current-phase relation stems from the fact that 
kx in eq. ([3|) becomes complex number in the presence of 
rrix. The position of the present CMMs in S/FI/S junc- 
tions is given by Ejj = ±Ej. The expression of E'^, can be 
explained by the hybridization of two CMM formed at the 
left S/FI interface and right FI/S interface. The formula 



FIG. 4: (Color online) Contour plot of CMM energy level Ej 
as a function of and cp for (a)mcc = 0, {h)mx = OArriz, 
{c)mx = —OArriz. The resulting Josephson current in S/FI/S 
junctions is plotted in (d) with a: mx/rriz = 1, b: mx/rriz — 

0. 4 and c: rrix/mz = —0.4. In all panels, we choose rrizd/vF = 

1, ii/rriz — 1 and rriy/mz — ^. T — 0.05Tc with transition 
temperature Tc- 

Differently from cr for the Andreev reflection, / is al- 
most independent of rriy because the contributions to the 
change in / almost cancel between the two CMMs with 
opposite chiralities. On the other hand, rrix seriously 
influences / as the effective vector potential which di- 
rectly enters into the phase of the Josephson coupling 
if. Here, the absence of the small factor of e/c, which 
reduces the coupling to the magnetic field, makes the 
tuning of the current-phase relation much easier. Re- 
flecting the dependence of Ej^ the resulting / becomes 
zero at neither (/p = nor ip = ±7r (curves h and c in 
Fig. 4(d)) due to the presence of the phase shift 25. Up 
to now, there have been many studies about Josephson 
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current in ferromagnet junctions [20], where the value 
of phase shift becomes neither nor tt in usual cases. 
The intermediate phase shift has been predicted in a few 
cases even in the absence of Major ana Fermion, i) spin- 
singlet 5- wave / spin-triplet superconductor junction [21], 
ii) even- frequency / odd- frequency superconductor junc- 
tion [22i], and iii) ferromagnet junction with spin-singlet 
s-wave superconductor in the presence of the spin-flip 
scattering or spin-orbit coupling 1^ . It is remark- 
able that this anomalous current-phase relation appears 
simply controlling one magnetization vector without us- 
ing unconventional pairing in the present model. It is 
also noted that the resulting current-phase relation can 
be continuously tuned by the change of rrix similar to 
the control of the magnetization vectors at the interfaces 
of ferromagnet junction [23]. We hope this anomalous 
current-phase relation will be detected experimentally in 
SQUID. 

In conclusion, we have studied the charge transport 
properties of normal metal (N) / ferromagnet insulator 



(FI) / superconductor (S) junction and S/FI/S junction 
formed on the surface of three-dimensional (3D) topo- 
logical insulator, where chiral Majorana mode (CMM) 
exists at FI/S interface. We have found that CMM gen- 
erated in N/FI/S and S/FI/S junctions are controlled by 
the direction of the magnetization m in FI region very 
sensitively. Since metallic surface state of 3D topological 
insulator has been observed experimentally, our proposed 
structure will be accessible in the near future [25|. Our 
results give a guide to innovate novel direction of super- 
conducting spintronics. 

Note added: After a submission of thie paper, a re- 
alization of the Majorana fermion in superconductor / 
semiconductor heterostructure has been predicted in the 
presence of interplay of spin-orbit coupling and exchange 
field [26]. 
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